Sub-Poissonian atom number fluctuations using light-assisted collisions 
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We investigate experimentally the number statistics of a mesoscopic ensemble of cold atoms in a 
microscopic dipole trap loaded from a magneto-optical trap, and find that the atom number fluc- 
tuations are reduced with respect to a Poisson distribution due to light-assisted two-body collisions. 
For numbers of atoms N > 2, we measure a reduction factor (Fano factor) of 0.72 ± 0.07, which 
differs from 1 by more than 4 standard deviations. We analyze this fact by a general stochastic 
model describing the competition between the loading of the trap from a reservoir of cold atoms 
and multi-atom losses, which leads to a master equation. Applied to our experimental regime, this 
model indicates an asymptotic value of 3/4 for the Fano factor at large A'^ and in steady state. We 
thus show that we have reached the ultimate level of reduction in number fluctuations in our system. 

PACS numbers: 05.40.-a,03.65.Ta,34.50.Rk,42.50.Lc 



There is a growing interest in the study of mesosco- 
pic systems containing between 10 — 100 particles. For 
example, mesoscopic ensembles of ultra-cold atoms could 
be a useful tool for quantum information processing (e.g. 
P, Q) or for precision measurements beyond the stan- 
dard quantum limit [3, '3| . They are also a test bed for 
the investigation of many-body correlated quantum sys- 
tems pI and for the study of collective effects such as 
super-radiance (e.g. (i'l). All these applications require 
the precise knowledge of the distribution of the number 
of atoms as the properties of these finite size samples 
are governed by their statistical nature. In particular 
the knowledge of the variance of the number distribu- 
tion is important. In this paper we show experimentally 
and theoretically that one of the conceptually simplest 
mesoscopic systems, namely a few cold atoms in a tight 
dipole trap being loaded from a cold atomic cloud [THQl, 
already exhibits non trivial sub-Poissonian statistics. 

The preparation of an atomic sample with non- 
Poissonian atom number distribution requires a non- 
linear mechanism usually provided by interactions bet- 
ween ultra-cold atoms. For example the dispersive s-wave 
interaction was used to reduce the relative atom num- 
ber fluctuations between the sites of a two-well poten- 
tial 0, 11 1 or of an optical lattice [13] ■ This led in par- 
ticular to the study of the Mott transition [H , 14 j . This 
interaction was also used to demonstrate reduced atom- 
number fluctuations in a single tight dipole trap [l5| . 
Recently, the production of a sample with definite atom 
numbers was demonstrated using the Pauli blockade [l6| . 

Inelastic collisions between ultra-cold atoms can also 
provide the non-linearity required to modify the atom 
number statistics, as shown recently [13, those ex- 

periments, three-body inelastic collisions induce losses in 
an initially trapped sample of 50 — 300 atoms at or close 
to quantum degeneracy, and the ever decreasing frac- 
tion of remaining atoms exhibits reduced number fluc- 
tuations with respect to a Poisson distribution. Here, we 



consider theoretically a different and yet more general re- 
gime where the trap continuously experiences the inter- 
play between a loading process from a reservoir of laser- 
cooled atoms and strong inelastic p-body losses {p > 1), 
and we investigate experimentally the case p = 2 where 
the losses are due to light-assisted collisions. This situa- 
tion is used elsewhere to, e.g., produce a single atom 
source [3, H^l- There, one operates in the "collisional 
blockade" regime where the loading rate is such that 
the microscopic trap contains one or zero atom with 
equal probabilities {{N) — 0.5) and the atom number 
distribution is maximally sub-Poissonian with variance 
AN^ = 0.5 (N) [iH. 

Here, we explore the regime where (N) goes beyond 
0.5 in steady state as we increase the loading rate. In 
practice, we prepare a thermal ensemble of up to 10 cold 
atoms at a temperature of ~ 100 fiK in a microscopic 
dipole trap. We observe that the atom number distribu- 
tion remains sub-Poissonian and that the reduction in 
number fluctuations with respect to the Poisson distri- 
bution, AN'^/{N), is locked to a constant value of 0.75 
for (N) > 2, a fact that was overlooked so far. To explain 
this fact, we use a microscopic approach that takes into 
account the stochastic nature of the competing loading 
and loss processes, and we calculate the atom number 
distribution at any time of the system evolution. We do 
so by solving a master equation both numerically and 
analytically, and find good agreement with the average 
result of a Monte Carlo approach where we study the 
individual behavior of atoms. Using this general theore- 
tical approach, we analyze our data and find that we have 
reached experimentally the ultimate level of reduction in 
atom number fluctuations that one can expect in a dipole 
trap operating in our regime {p — 2). The formalism pre- 
sented in this paper is applicable to any system where 
a random loading process competes with a p-body loss 
process, whatever its nature and whatever the number of 
atoms. 
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To study the number statistics of a few atoms in the 
presence of competing random processes, we implemen- 
ted the fohowing experiment (details can be found in 
Ref. (i). First, we produced a microscopic optical dipole 
trap at 850 nm by sharply focusing a laser beam [22| . 
We then loaded this trap from a magneto-optical trap 
(MOT) of *^Rb atoms surrounding the region of the di- 
pole trap. Atoms enter the dipole trap randomly, are 
trapped thanks to the cooling effect of the MOT beams, 
and are expelled from the trap due mainly to inelastic 
two-body collisions assisted by the near-resonant light 
of the cooling beams and, in a minor extent, to colli- 
sions with the residual background gas in the chamber 
(one-body losses). We measured elsewhere the two-body 
and one-body loss constants to be /?' ~ 500 (at.s)~^ and 
7 ~ 0.2 respectively [2^. The actual value of the 
loading rate R is proportional to the MOT local density 
in the vicinity of the microscopic trap, which is the pa- 
rameter that we vary. For values of i? 3> /?' the mean 
number of trapped atoms in steady state exceeds unity 
and is (A^)st = while for i? <C 7 it goes to zero 

as {N)st = R/j- The intermediate regime corresponds 
to the "collisional blockade" regime where (A^)st =0.5. 
Experimentally we operate at (A^)st ^ 2 in the following. 

To get information on the number distribution of 
atoms in the dipole trap in steady state, we release the 
atoms from the trap and probe them with a pulse of reso- 
nant light. Using an intensifier to amplify single photon 
events above the noise of our CCD camera^e count the 
detected fluorescence photons individually [9| . This num- 
ber is proportional, on average, to the number of atoms N 
in the trap before release. We build up the number distri- 
bution of counted photons by repeating this loading and 
probing experiment about 100 to 1000 times. Knowing 
the response of our imaging system to one atom exactly, 
we extract from the photon distribution the mean {N)st 
and the variance AiV^ of the atom number distribution 
in steady state, and calculate the corresponding Fano fac- 
tor F = AN^/{N)st- The data shown in Fig. [T] indicate 
a clear reduction of the atom number fluctuations with 
respect to a Poisson distribution for {N)st > 2 with a 
mean F = 0.72 and a total uncertainty (1 standard de- 
viation) of 0.07. This uncertainty is obtained by adding 
quadratically the statistical (type A) uncertainty of 0.05 
(deduced from the rms dispersion of the data) and the 
systematic (type B) uncertainty of 0.04, which we esti- 
mated in previous work 9]. 

Qualitatively, this reduction can be understood as fol- 
lows. If the losses were governed by random one-body 
events, e.g. background gas collisions, the trade-off bet- 
ween the random loading of the trap and the losses would 
result into a Poisson distribution with mean atom num- 
ber (TV) St in steady state. If the losses now involve higher- 
body processes {p > 2) the loss rate varies as the number 
of p-uplets in the A'^-atom ensemble, i.e. increases non- 
linear ly with TV. For a given mean atom number this 




Figure 1: Fano factor F versus the average number of atoms 
(A'^)st in the microscopic dipole trap in steady state. Square 
and triangles : experimental data collected in the "collisional 
blockade" regime {{N)st = 0.5) and beyond ((A'')st > 0.5). 
Solid line : model based on a stochastic process (see text) with 
= 500 (at.s)~^ and 7 — 0.2 s~^. Dashed line : theoretical 
limit F — 3/4: for {A^)st 3> 1. f : rms dispersion of the data 
collected beyond the collisional blockade regime. 



leads to the number distribution being narrower than a 
Poisson distribution, as the losses are more efficient on 
the high-A side of the distribution. 

To explain quantitatively the sub-Poissonian behavior 
of the atom number distribution, we use the following sto- 
chastic model that takes into account the three random 
processes involved, i.e the loading, the two-body losses, 
and the one-body losses. We consider the evolution in 
time of the probability PAr(i) to have A^ atoms in the 
dipole trap. To calculate the probability p^it + dt), we 
sum the contributions of all channels associated to the 
random processes mentioned above that lead to having 
A atoms in the trap ai t + dt, given that the trap could 
possibly be filled with either A^ - 1, A^, A -I- 1, or A^ -h 2 
atoms at time t. The probability that a loading event 
occurs in the time interval dt when there are already A 
atoms in the trap is Rdtpf^{t). Similarly, the probability 
that a loss event occurs during dt is jN dtp^it) for one- 
body events, and /3' dtpN(t) for two-body events. 
We obtain eventually the following equation : 

N(N - 1) 

PN{t + dt)= pN{t) (1- [i? + 7iV + /?'^^ '-]dt) 

+PN-i{t) Rdt 
+PN+l{t) l{N + l)dt 

(iV + 2)(iV + l) 
+PN+2[t) 13 dt. (1) 

Taking the limit — > 0, Eq. ((T|) yields the following 
master equation that rules the evolution of pn (t) in time : 

^^"^ R (E-i-l)M 

-^7 {I]~1)[Npn] 



dt 



+ f3' (E2-1)[ 



N{N - 1) 



-Pn\ 



(2) 



where E is the "step operator" defined by its effect on an 
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Figure 2: Set of numerical solutions {p]v(t)} of eqs. la- 
beled by N. Parameters : R = 6000 s^^ 7 = 0.2 s~^ and 
/3' = 500 (at.s)-i lead to {N)^t = 3.6 ~ ^/BJW- Inset : filled 
bars : the set of solutions in steady state {pjv(t = 3 ms)} (filled 
bars) is clearly sub-Poisson ; sticks : Poisson distribution with 
same mean value. 



arbitrary function f{N) : 

E[/(7V)] -/(7V + 1), E-i[/(A^)] (3) 

and 1 is the identity operator. Using eqs. 1^ we obtain 

the equation of evolution of the mean number of atoms 

00 

{N) = E : 

N=0 



d{N) 
dt 



R-j{N}- I3'{N}{{N} -1)~I3' AN^ . (4) 



When AA^^ = 0, we recover the phenomenological equa- 
tion sometimes used to describe the loading of a trap 
containing a small number of atoms [llll^], i.e. dN/dt = 
i?-7 N-P' N{N-l). When AA^^ ^ assuming 
a Poisson distribution) , eq. ^ also yields the widely used 
equation d{N)/dt = i? - 7 (iV) - /?' {N)'^ . However, wi- 
thout any a priori relation between AiV^ and (iV) , eq. (|4]) 
cannot be solved analytically. 

To calculate the first moments of the number distri- 
bution, we used three different approaches. First, we sol- 
ved numerically eqs.© using the boundary conditions 
PAr(O) = SNfi and, for N » (7V)st, pjv(i) = 0. As an 
example. Fig. [2] illustrates the time evolution of the pro- 
babilities PN{t) for parameters leading to {N)st = 3.6. 
The number distribution is found to be sub-Poissonian, 
with F = 0.74. By varying the loading rate the same ap- 
proach yields the distribution for any value of {N)st- We 
analyze the case where 7 <C /?' in the following. When 
(TV) St ^ 0.5 we find, as expected, that the distribution is 
close to a Poisson law as one-body losses then dominate 
two-body losses (see Fig. |31^a)) [27[. The presence of two- 
body processes induced losses of atom pairs leads to a 
sub-Poissonian behavior that is maximal for (iV)st =0.5, 
corresponding to Po = Pi = 0.5. While this regime has 
been described before (see Ref. (HI), the numerical ap- 
proach predicts that atom number fluctuations do not be- 
come Poissonian for larger numbers of atoms. In fact, the 
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Figure 3: Theoretical results obtained by different ap- 
proaches. Solid line : numerical solution of eqs. ([2]) for our 
experimental parameters (7 — 0.2 s~^, /3' ~ 500 (at.s)~^) ; 
squares : Monte-Carlo simulation, (a) Dependence of the Fano 
factor on {N}af (b) Example of calculated atom number 
distribution in steady state, with R = 5 x 10^ s~^, yiel- 
ding {N)st = 32. The numerical solution is indistinguishable 
from the Gaussian solution of eq. ((SJ. We find a Fano fac- 
tor F = 3/4. Dotted line : Poisson distribution with the 
same mean value. Dashed line in (a) : numerical solution 
of eqs. ([2]) when two-body processes induce the loss of one 
atom only from the trap (parameters : 7 = 5 x 10^'^ s~^, 
/3' ~ 500' (at.s)-^). In the limit 7/^' 0, (iV>st can be lo- 
cked to 1 in a fully deterministic way (AA'^'^ — 0). 



Fano factor reaches an asymptotic value of 0.75 as soon 
as (A^)st ^ 2, corresponding to a reduction of —1.25 dB 
with respect to the Poisson case. The numerical predic- 
tion reproduces well our data, as shown in Fig.[T] 

The second approach to solve eqs. ^ is analytical. It is 
valid for (N)st ^ 1 only and follows closely the approach 
of Ref. [25| . We first re- write the master equation into a 
dimensionless rate equation : 



dp 



N 



dT 



t)[p 



-Pn], 



(5) 



where the one-body loss term of eq. ^ has been neglec- 
ted (following i? ^ 3> 7) and r = t^JR(3' is a dimen- 
sionless time variable. Since the number distribution in 
steady state is expected to be peaked around (A^)st with 
a width on the order of {N)st, we consider the number 
of trapped atoms at time t as a stochastic quantity of 
the form 



7V(r) = (iV),t 0(r) + /(A% ^(r) 



(6) 



where ^(r) is a stochastic variable with mean (^) (r) = 
and an amplitude of '-^ 1, and 0(t), also on the order of 



4 



1 , is a deterministic and slowly varying function of time 

(0(r) = (N) (t) / (N) St) ■ We then consider the probability 

P(^,t) — pn{t) that N atoms are in the trap at time 

T. Since PN+k{r) = P{C + ,t) and (A^)st > 1, we 

V \")st 

replace pAr+fe(T) in eq.(l5]) by aTaylor expansion of P(^, r) 
in powers of 1/a/ (iV)st- Replacing pn{t) by P{^,t), the 
time derivative dpN/dr becomes drP — VNst </* d^P [28| . 
Identification of the power terms in the expanded master 
equation then yields the following equations that rule the 
evolution of (p and P in time : 

^ = (7) 

drP = 2^ d^i^P) + i(2^2 _^ ^) Q2p 

Eq. ([8]) is a linear Fokker-Planck equation with time de- 
pendent coefficients, the steady state solution of which 
is Gaussian [2^. Finally, using eq. dS]) we find that (^^) 
evolves in time according to : 

^ = -40(a + (l + 20^). (9) 

Equations (O and ^ allow us to calculate the evolution 
of the Fano factor in time, F{t) — i^"^) {t) / 4>{t) . In parti- 
cular, in steady state, (f) = 1 and F = 3/4. This analytical 
finding is in excellent agreement with our numerical so- 
lution (see Fig. |3]). Besides, we find that the analytical 
result is valid for atom numbers as small as ^ 2. 

Finally, we cross-checked our theoretical results by 
a Monte-Carlo simulation, where we calculate at each 
time increment the survival probabilities of individual 
atoms to the various random events involved in the pro- 
blem [2^. By averaging over many atomic histories, we 
reconstructed atom number distributions and found Fano 
factors in very good agreement with those presented 
above (see Fig. ([3])), which validates the master equation 
approach. 

In conclusion, we discuss our experimental findings and 
theoretical approaches from a more general perspective. 
First, the observed reduction in number fiuctuations is 
due to loss terms that vary non-linear ly as iV, and is thus 
intrinsically robust to losses (provided p and the loading 
rate remain constant). In our case, we reached experimen- 
tally the ultimate level of reduction (—1.25 dB) predicted 
by theory when a loading mechanism competes with a 
two-body non-linearity leading to the loss of atom pairs, 
no matter the underlying mechanism (light-assisted col- 
lisions, hyperfine changing collisions...). More generally, 
the exact level of reduction achievable depends on p, on 
the number of atoms being lost after a p-body process, 
and on the presence (or the absence) of a loading mecha- 
nism. When R ^ 0, the analytical approach explained 
above can be generalized and yields a Gaussian atom 
number distribution in steady state (/'(C, r) evolves ac- 
cording to a Fokker-Planck equation similar to eq. ([5])). 
For p-body processes leading to losses of p-uplets, one 



finds an equation similar to eq. ([9]) and F = 5(1 + ^) 
in steady state. When R = 0, slightly better levels of re- 
duction can be achieved, as F — p/ {2p — 1). This was re- 
cently demonstrated in the case of three-body losses [3] ■ 
Finally, we extended our approach to the case where two- 
body collisions lead to the loss of one atom only from 
the trap [2^ , as is the case for elastic collisions induced 
evaporative losses and for some light-assisted loss me- 
chanisms. Such mechanisms have been used recently to 
produce near-deterministically a sing le atom source for 
quantum information processing [26|. Taking these me- 
chanisms into account, our theoretical approach predicts 
that fiuctuations fully vanish when only 1 atom is left in 
the trap in the absence of one-body decay (see Fig.[3l[a)), 
i.e. that a robust and fully deterministic preparation of 
single trapped atoms is in principle possible. 
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